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I 1 INTRODUCTION 



ABSTRACT 

Weak gravitational lensing provides a direct way to study the mass distribution of 
clusters of galaxies at large radii. Unfortunately, large scale structure along the line of 
sight also contributes to the lensing signal, and consequently affects the measurements. 
We quantify the effect of distant uncorrected large scale structure on the inferred 
mass profile of clusters as measured from weak lensing. We consider NFW profiles, 
and find that large scale structure is a major source of uncertainty for most practical 
situations, when a model, with the mass M200 and the concentration parameter c 
as free parameters, is fit to the observations. We find that the best constraints are 
found for clusters at intermediate redshifts (z s» 0.3). For a cluster at z = 0.3, optimal 
results are obtained when the lensing signal is measured out to 10 — 15 arcminutes. 
Measurements at larger radii do not improve the accuracy with which the profile can 
be determined, contrary to what is expected when the contribution from large scale 
structure is ignored. The true uncertainties in M200 and the concentration parameter 
c are ~ 2 times larger than when distant large scale structure is not included in the 
error budget. 
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Dark matter plays an important role in the formation of 
structure in the universe. However, little is known about 
the nature of dark matter, but it has become clear that 
the simplest paradigm, that of collisionless cold dark matter 
(CDM), can explain the observed structure in the universe 
quite well. Numerical simulations, which provide a powerful 
way to study the formation of structure in the universe, 
indicate that CDM gives rise to a particular density profile 
(e.g., Dubinski & Carlberg 1991; Navarro, Frenk, & White 
1995; Navarro, Frenk, & White 1997, hereafter NFW; Moore 
et al. 1999). In these simulations, the NFW profile appears 
to be an excellent description of the radial mass distribution 
for halos with a wide range in mass. 

The density profile on small scales, however, remains 
controversial as other groups find different results (e.g., 
Ghigna et al. 2000). In addition, realistic simulations should 
include the effect of baryons, which complicate matters even 
further. On the other hand, there is good agreement for the 
density profile on large scales, and a comparison of the pro- 
files of real objects with the predictions provides an impor- 
tant test of the assumption that structures form through 
dissipationless collapse. 



Clusters of galaxies are particularly well suited for such 
a comparison, because their formation can be simulated 
fairly well. Unfortunately, because clusters are dark matter 
dominated, it is difficult to measure the mass distribution 
at large radii from the cluster centre. Dynamical estimators 
can be used, but assumptions have to be made about the 
shape of the cluster, its dynamical state, and the velocity 
anisotropies (e.g., van der Marel et al. 2000). 

A promising approach is that of weak gravitational lens- 
ing. The tidal gravitational field of the cluster distorts the 
images of distant background galaxies. This distortion can 
be measured and can be used to reconstruct the projected 
mass density of the cluster (e.g., Kaiser & Squires 1993). 
The advantage of weak lensing over dynamical methods is 
that the results do not rely on the dynamical state or the na- 
ture of the deflecting matter. More importantly, it does not 
require a tracer of the gravitational potential (which limits 
many dynamical methods) , and as a result the lensing signal 
can be measured out to large radii from the cluster centre. 

The weak lensing signal induced by galaxy clusters has 
been measured for a large number of systems (e.g., Bon- 
net, Mellier & Fort 1994; Fahlman et al. 1994; Clowe et 
al. 2000; Hoekstra et al. 2002a; Squires et al. 1996). These 
early studies were limited by the field of view of the CCD 
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cameras used in the observations. With the recent introduc- 
tion of wide field imagers it is now possible to measure the 
lensing signal out to large radii (Clowe & Schneider 200f, 
2002) or study the mass distribution in superclusters (e.g., 
Kaiser et al. f998; Gray et al. 2002). Many other weak lens- 
ing applications are discussed in Hoekstra, Yee & Gladders 
(2002b). 

The lensing signal at large radii is low, and special care 
has to be taken to correct for observational distortions, such 
as PSF anisotropy, seeing and camera shear. These system- 
atic signals typically have amplitudes comparable to the 
lensing signal one is interested in. Fortunately, detailed stud- 
ies have shown that the lensing signal can be recovered with 
great accuracy (e.g., Hoekstra et al. f998; Erben et al. 200f ; 
Bacon ct al. 200f). 

However, another source of error is usually ignored: the 
contribution of large scale structure. The weak lensing sig- 
nal is sensitive to all matter along the line of sight. In gen- 
eral, cluster studies assume that the signal introduced by 
the cluster is so dominant that the contributions from other 
structures can be neglected. 

The effect of local large scale structure (such as the fila- 
ments connecting clusters and groups) on cluster mass esti- 
mates has been studied through numerical simulations (Cen 
1997; Metzler et al. 1999; White et al. 2002). These studies 
show that local structures introduce both a bias and addi- 
tional noise in the mass measurement. Estimates of cluster 
mass-to-light ratios are affected less by these structures, as 
both the mass and the light are expected to probe the same 
structure. 

The study of the contribution from distant large scale 
structure (not correlated with the cluster) requires large cos- 
mological simulations, as one needs to account for all mat- 
ter along the line of sight. Hoekstra (2001) demonstrated 
how the contribution of distant large scale structure can 
be computed analytically. Hoekstra (2001) showed that the 
large scale structure introduces noise in the cluster mass 
measurement, but does not bias the result. The results pre- 
sented by White et al. (2002), based on detailed numerical 
simulations, are in good agreement with the findings from 
Hoekstra (2001). 

It is important to note that the distant large scale struc- 
ture does affect the determination of the cluster mass-to- 
light ratio. Furthermore, the distant large scale structure 
provides a fundamental limit to the accuracy of weak lensing 
mass estimates: deeper observations reduce the statistical 
uncertainty (because of the intrinsic shapes of the sources) , 
but increase the noise from large scale structure (Hoekstra 
2001). 

These studies have demonstrated that the large scale 
structure can be an important source of noise for the mass 
determination of individual objects. The cluster mass is an 
integral over the cluster profile, and therefore large scale 
structure will have an impact on the determination of cluster 
mass profiles. 

Local large scale structure and substructure in the clus- 
ter will change the profile. This was studied by King, Schnei- 
der & Springel (2001) using a numerical simulation of a clus- 
ter. They found that substructure increases the error on the 
parameters by only ~ 3%. Doubling the amount of substruc- 
ture resulting in ~ 10% larger errors. Based on these results, 



King et al. (2001) concluded that weak lensing studies of 
clusters can be of great use to study mass profiles. 

However, King et al. (2001) did not consider the contri- 
bution of distant large scale structure to the error budget. 
In this paper, we use the approach introduced by Hoekstra 
(2001) to quantify the contribution of distant large scale 
structure on measurements of cluster density profiles. 

The structure of the paper is as follows. In Section 2 we 
calculate the contribution of distant large scale structure to 
the lensing signal. The lensing properties of the NFW profile 
are discussed in Section 3. The effects of distant large scale 
structure are studied in Section 4, where we consider clusters 
with different masses and redshifts. 



2 NOISE FROM LARGE SCALE STRUCTURE 

In the weak lensing regime, the coherent distortion of the im- 
ages of the background galaxies provides a direct measure 
of the lensing shear 7 as a function of position. The pro- 
jected mass distribution of a cluster can be reconstructed 
in a parameter- free way from the observed shear field (e.g., 
Kaiser & Squires 1993). The resulting mass map can be used 
to derive the cluster mass profile. However, the mass recon- 
struction is a non-local operation, and as a result the noise 
in the map varies as a function of position. 

Although one can estimate the noise in non-parametric 
reconstructions from bootstrap simulations, it is more con- 
venient to use parameterised models which are fitted to the 
observations. In the latter case the error properties are well 
understood. To study the distribution of matter in the clus- 
ter, 2-dimensional models can be fitted. Another approach 
is to consider the azimuthally averaged tangential shear as a 
function of radius from the centre. The latter is particularly 
useful to estimate the mass of the cluster (e.g., Hoekstra et 
al. 1998) or study the mass profile (e.g., Clowe et al. 2000; 
Clowe & Schneider 2001; Clowe & Schneider 2002; Hoekstra 
et al. 2002a). 

Hoekstra (2001) showed how distant large scale struc- 
ture affects the mass estimate when a singular isothermal 
sphere (SIS) model is fitted to the tangential shear profile. 
In this paper we apply the approach presented in Hoekstra 
(2001) to estimate how large scale structure affects the mea- 
sured tangential shear profile. This allows us to efficiently 
simulate observations, and quantify the effect of distant large 
scale structure on the derived cluster profile. We focus on 
the NFW profile, which is discussed in Section 3, although 
the procedure can be readily applied to any model. 

The quantity of interest is the tangential shear aver- 
aged in an annulus ranging from 61 to 62: (jt) {81, 82)- We 
will show that the tangential shear in an annulus is a fil- 
tered measurement of the surface mass density, i.e., it is a 
particular choice of the aperture mass statistic M ap (e.g., 
Schneider et al. 1998). The aperture mass is defined as 

M ap (#)= / d 2 m(W)K(-ff), (i) 

Jo 

where U(8) is the weight or filter function. Provided that 
U(8) is compensated, i.e. 

[ dd tf?7(tf) = 0, (2) 
Jo 



© 0000 RAS, MNRAS 000, 000-000 



Mass profiles from weak lensing 3 



the aperture statistic can be related to the tangential shear 
7t. Unfortunately 7t cannot be observed directly, but one 
can measure the distortion gr = 7t/(1 — «)• In the weak 
lensing regime, where k< 1, one can use 7t ~ gr, which is 
what we do throughout this paper. Hence, 



M ap (0) 



Jo 



Q(0) is related to U(Q) through 



Q(0) 



2 
02 



dddUiti) - U{6). 



(3) 



(4) 



Hence, similar to the mass estimate of the cluster 
(Hoekstra 2001), the tangential shear in an annulus can be 
written as an aperture mass, with the filter function Q{9) 



Q(\e\) 







{0i < \e\ < ft,) 

(elsewhere) 



and the corresponding f/(|0|) is 



u(\e\) 




(9 < 0i) 

(8i < W < 
(9 > 9 2 ) 



The coefficients on are defined as 
2[ln(6> 2 )-ln(6> 1 )] 



ai = 



Ct2 = 



<9l 
-2 



and 



a?, 



21n(0 2 ) - 1 



(5) 



(6) 



(7) 



(8) 



(9) 



<9l-9\) ■ 

Under the assumption that the cluster is the only lens- 
ing structure in the field, the measured (7t)(0i,02) as a 
function of the distance to the cluster centre can be readily 
used to study the cluster mass profile: the parameters de- 
scribing the cluster can be inferred from a least squares fit 
of the model to the data. However, the observed weak lens- 
ing signal is sensitive to all matter along the line of sight, 
and other structures (the distant large scale structure) will 
introduce an additional signal, and the observed signal is 
the sum of the contribution from the cluster and large scale 
structure 



(7T)obs = (7T)clustor + {~/t)lSS- 



(10) 



The aperture mass measures a density contrast, and 
therefore the expectation value for the contribution by un- 
corrected large scale structure vanishes (i.e., (Mlss) = 0). 
Thus on average, distant structures do not bias the mea- 
surements, but they introduce and additional statistical un- 
certainty in the measurement of size cr 2 ss = ((7t , )lss) 1//2 - 
However, because of non-linear gravitational evolution, the 
distribution of M ap will be skewed (e.g., Bernardeau et al. 
1997; Schneider 1998). In the analysis we will ignore the 
skewness, and assume that M ap follows a normal distribu- 
tion. 

The uncertainty in the measurement of the average tan- 
gential shear a ohB in a bin ranging from 0i to 92 is 



2 2 2 

Cobs = Ogal + 0- LSS , 



(11) 




Figure 1. The function g 2 (l, 6i, #2) as a function of I, for 0\ = 15' 
and 82 = 15'.5. The filter function declines rather slowly (oc 1/Z). 



where <r| al corresponds to the uncertainty introduced by the 
intrinsic shapes of the galaxies and the noise in the shape 
measurement. It depends on the number density of sources 
n and the scatter in the ellipticities cr e 



2 



<9l 



9?)n 



(12) 



For the dispersion in the shapes of the sources we adopt a 
value <r e = 0.25, similar to what is observed in deep HST 
images (e.g., Hoekstra et al. 2000). Depending on the image 
quality and depth of the observations, the dispersion is likely 
to be larger for faint galaxies in ground based data, because 
of the correction for the circularization by the PSF. The 
number density of galaxies with 20 < R < 26 is ~ 60 galax- 
ies arcmin -2 in space based observations. In ground based 
data, however, the effective number density of galaxies is 
lower, because many of the faint galaxies have sizes compa- 
rable to the PSF, and are not useable in the weak lensing 
analysis (i.e., these galaxies have cr £ 3> 0.25). We therefore 
adopt a typical value of 30 galaxies arcmin -2 (e.g., Bacon 
et al. (2002)). 

The variance in the aperture mass caused by large scale 
structure is given by (Schneider et al. 1998) 



&lss(9i,9 2 ) 



= 2tt I 
Jo 



dl lP K (l)g 2 (l, 0i,0 2 ), 



(13) 



where the effective projected power spectrum P K (l) is given 
by: 



I 



w (14) 



Here w is the radial coordinate, a(w) the cosmic scale factor, 
and /k(w) the comoving angular diameter distance. W(w) 
is the source averaged ratio of angular diameter distances 
Dis/Ds for a redshift distribution of sources p m (w): 
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Figure 2. The dispersion oxssC?) m the averaged tangential 
shear introduced by distant large scale structure. The solid lines 
shows the results for sources with 20 < R < 26 and the dashed 
line corresponds to sources with 20 < R < 24. The effect is smaller 
for brighter (lower redshift) sources, but in this case the cluster 
lensing signal is lower as well. 



Figure 3. The cross-correlation coefficient r for 8 = 15' (solid 
line), 8 = 10' (dotted line) and 8 = 5' (dashed line). The mea- 
surements on at small radii are not strongly correlated, but at 
large radii the correlations between bins is rather strong. 



T J>/ x f WH , i , tsfK{w'-w) 

W(w) = I dw p w {w ) 



(15) 



f K (w>) ■ 

The function g(l, 61,62) in equation [l^ depends on the filter 
function U(\6\) as 



9(1, 01,02) 



d4>4>U{4>)Mi4>). 



(16) 



As has been demonstrated in Jain & Seljak (1997) and 
Schneider et al. (1998) it is necessary to use the non- linear 
power spectrum in equation 14. This power spectrum can be 
derived from the linearly evolved cosmological power spec- 
trum following the prescriptions of Hamilton et al. (1991), 
and Peacock & Dodds (1996). 

For our choice of filter function C/(|#|) we obtain 



9(1, 



n,02 = 



21n(6»i)\ 0iJi(Z6>i) 



<oi - el) 



+ a3 6 2 Jl(l6 2 ) 1 ,/ 1 „.ll l (..)./n(/o|.(l7l 



Figure [l] shows g 2 (l, 61, 62) as a function of / for 61 = 15' 
and 62 = 15'.5. These results show that g 2 declines rather 
slowly (oc 1/1) with increasing I, and therefore probes a large 
range in I. 

For the analysis in this paper it is convenient to consider 
the aperture mass statistic 5"lss(#); which corresponds to 
the the average tangential shear measured in an annulus 
[6 - 86/2, 6 + 86/2]. Provided that 86/6 < 1, one obtains a 
simple form 

Mid) 



g(l,e,dd)^g(ie) = 



271 



(18) 



which is independent of 86, and is an excellent approxima- 
tion to the exact calculations. In the remainder of the paper 
we will adopt a value of 86 = 0'.5. 

Figure | shows the dispersion 5"lss(6 i ) in the averaged 
tangential shear introduced by distant large scale structure. 
To derive these results, we have used the redshift distribu- 
tion presented in Hoekstra (2001), which are based on pho- 
tometric redshift studies of the Hubble deep fields. The solid 
lines shows the results for sources with 20 < R < 26 and the 
dashed line corresponds to sources with 20 < R < 24. 

The number density and mean source redshifts are 
higher in deeper observations. Consequently, the lensing sig- 
nal is higher, and the statistical noise is lower. However, as 
Figure H shows, the effect of large scale structure is also 
increased, and the actual accuracy of the measurements im- 
proves only marginally (also see Hoekstra 2001). 

To study the effect of distant large scale structure on the 
measurement of the cluster mass profile, we create simulated 
tangential shear profiles. We include the noise caused by the 
intrinsic shapes of the background galaxies, and the noise 
introduced by intervening large scale structure. 

The statistical error a ga i in each bin, caused by the 
intrinsic shapes of the sources, is uncorrelated with the 
other bins. The noise from large scale structure at various 
scales is correlated. To properly account for this, we need 
to compute not only olss but also the cross correlation co- 
efficient between bins. The cross-correlation between bins 
(5"lss(0)5'lss(6 i — "&)) is calculated using 



<ft.Bs(0i)ffL8s(ty)> = 2vr / dllP K (l)g(l,6)g(l,6-tf). (19) 
Jo 

and the cross-correlation coefficient r(6; 6 — $) is 
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r(8;8--d) 



(^LSs(flKss(fl-#)) 



(20) 



The cross-correlation coefficients r for three different 
values of 8 are shown in Figure ^| The measurements on 
small scales have small range correlations, whereas at large 
radii the contribution of large scale structure introduces 
large scale correlations. 

Correlated random noise c(8) can be obtained by 
smoothing of the white noise n(8) with a an appropriate 
smoothing kernel K(8 — •&) 



c(8) 



K(8 - 0)n(0)d#, 



(21) 



where n(9), c(8) are random functions with correlators 
{n(6)n(0')) =6(0-8'), (22) 



and 



(c(8)c(8')) = / K (8 - $)K(8' - ■&)<!■&. 



(23) 



For n(8) we take a random number with a Gaussian 
distribution with a dispersion of unity. The smoothing kernel 
is readily obtained from the cross correlation coefficient. The 
correlation coefficient is the convolution of the smoothing 
kernel with itself. The smoothing kernel is normalized such 



that<c 2 (0)> = a 2 ss (0). 



3 NFW PROFILE 

Numerical simulations have indicated that dark matter ha- 
los originating from dissipationless collapse of density fluc- 
tuations can be described by a universal density profile 
(e.g., Dubinski & Carlberg; Navarro, Frenk, & White 1995; 
Navarro, Frenk, & White 1997). The NFW density profile is 
given by 

S c p c 



p(r) 



(r/r s )(l+r/r s ) 2 ' 



(24) 



where p c is the critical density of the universe at the redshift 
of the halo. The overdensity of the halo is parameterized by 
5c, which is related to the concentration parameter c through 

* e = gg° ^ (25) 



3 ln(l + c) - c/(l + c) ■ 

The scale radius r s is the characteric radius of the halo, 
which depends on the virial radius r2oo and the concentra- 
tion parameter c as r s = r2oo / c. The virial radius is defined 
as the radius where the mass density of the halo is equal to 
200p c , and the corresponding mass M200 inside this radius 
is given by 

800tt 



M200 



-per 2oo- 



(26) 



Given the cosmology, redshift, and mass of the halo, 
r2oo follows immediately. The simulations have also shown 
that c and M200 are correlated, albeit with some scatter. 
Hence, given M200, the values of 5 C , and c can be computed 
using the routine CHARDEN made available by Julio NavarroQ 



The 



CHARDEN 



in be 

tittp : //pinot . phys .uvic . ca/ ~jf n/ charden 



obtained 



from 



In this case the halo profile is described by a single param- 
eter, e.g., the mass. However, one can also consider c as a 
free parameter in the model. Consequently comparison of 
the inferred value with the predictions of the numerical sim- 
ulations provides a direct observational test of our current 
understanding of structure formation. 

For the NFW profile, the tangential shear 7t as a func- 
tion of radius 8 is given by (Bartelmann 1996; Wright & 
Brainerd 2000) 



r s S c p c 



9<(x) 



,x <l 



1T (8) = { r^[«+41n(I)] ,x = l , 
r " s fs g> (x) ,x>\ 



(27) 



where x = 8/r s . The critical surface density E c is given by 
c 2 D 

Ec = feGAt' (28) 

where Di, D s , and Di s are, respectively, the angular diame- 
ter distances between the observer and the lens, the observer 
and the source, and the lens and the source. The functions 
g<(x) and g>(x) are defined as 



9<{x) = 



8arctanru / 

h -it hi 

a; VI - z 2 x 

4arctanh 1 



(x 2 - 1) 



(x 2 - 1) (1 - X 2 ) 



2\l/2' 



(29) 



and 



8 arctan . 



g> (x) 



X 2 \fl 



4 arctan . 



+ 



■In 



(1) 



(x 2 - 1) 



(x 2 



\3/2 



(30) 



Figure |4j shows the tangential shear for a cluster with an 
NFW profile. The cluster has a mass M200 = 1O 15 /i _1 M , 
and is at a redshift z = 0.3. The amplitude of the lensing 
signal corresponds to that for source galaxies with apparent 
magnitudes 20 < R < 26. The dashed line in Figure ^corre- 
sponds to a Singular Isothermal Sphere (SIS) model, which 
was matched to the NFW model at intermediate radii. 

In principle one can distinguish between the SIS model 
and the NFW profile by examining the lensing signal on ei- 
ther small or large scales. In practice, however, substructure 
in the cluster core, and contamination of the sample of back- 
ground galaxies by cluster members, complicates the small 
scale comparison. Furthermore, the profile on small scales 
is still debated (e.g., Ghigna et al. 2000). The lensing sig- 
nal at large radii can be measured accurately, provided that 
the shapes have been accurately corrected for observational 
distortions. 



4 CONSTRAINING MASS PROFILES 

To study cluster mass profiles, it is common to compare a 
parametrized mass model to the data. Ideally one would like 
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0.1 1 10 100 

/9[arcmin] 

Figure 4. Comparison of the tangential shear for a cluster with a 
NFW profile (solid line), and a cluster with a SIS profile (dashed 
line). The NFW cluster has a mass A/200 = 10 15 /i _1 A/ Q , and is 
at a redshift z = 0.3. The amplitude of the signal corresponds to 
that of source galaxies with apparent magnitudes 20 < R < 26. 
The amplitude of the SIS profile was matched to the NFW profile 
at intermediate scales. The arrow indicates a radius of lh" 1 Mpc. 



to minimize the number of parameters describing the model. 
In this section we consider the NFW profile (Eqn. 24), with 
M200 and c as the only two free parameters. In addition to 
these two parameters, one could introduce the slopes of the 
profile on small and large scales. 

To examine how well we can constrain the values of 
M200 and c, we create simulated tangential shear profiles 
as described in Section 2. The profiles include the statis- 
tical noise due to the intrinsic shapes of galaxies and the 
contribution from distant large scale structure. A typical re- 
alisation is presented in Figure ^a. The points indicate the 
simulated measurements, and the error bars reflect only the 
statistical error (the intrinsic shapes of the sources). The in- 
put NFW signal is indicated by the dashed line. The actual 
lensing signal, which is modified because of the large scale 
structure, isgiven by the solid line 

Figure |5p shows the difference between the "observed" 
signal and the input NFW profile. The large scale structure 
introduces significant deviations, which, in this particular 
case, systematically lower the signal on scales larger than 
10 arcminutes. The input profile has M200 = lO 15 ft -1 M0 
and c = 4.6. The best fit NFW profile to the profile 
with the contribution from large scale structure only (i.e., 
no statistical noise from the shapes of the sources) yields 
M200 = 8.7 x 10 14 A'/ Q and c = 5.1 when fitted out to 15 
arcminutes. When the profile is fitted out to 30 arcminutes, 
we obtain M200 = 7.4 x 10 14 A/ e and c = 5.5. Hence, the 
large scale structure alters the inferred profile of an individ- 
ual cluster, and the inferred values depend on the range in 
radius included in the fit. Note, however, that the ensemble 





M a0 o=i° 16h - 1 M G (a) . 




c = 4.64 


k 


z = 0.3 


- 
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10 20 30 
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Figure 5. (a) Realisation of the observed tangential shear as a 
function of radius from the cluster centre for a cluster at a redshift 
z = 0.3 and mass M200 = 10 15 /i — 1 Mq. The points include the 
noise caused by the intrinsic shapes of the background galaxies 
and the contribution from large scale structure. The error bars on 
the points indicate the (uncorrelated) statistical error in the mea- 
surements due to the shapes of the sources (adopting <r e = 0.25 
and ft = 30 arcmin - 2 . The model NFW profile is indicated by 
the dashed line. The cluster profile with the contribution from 
large scale structure added to the signal is indicated by the solid 
line, (b) The difference between the observed signal and the in- 
put NFW model. The large scale structure introduces significant 
deviations from the model. 

average value is not changed when a large sample of clusters 
is used. 

We now quantify the contribution of distant large scale 
structure to the total error budget. We start by consid- 
ering clusters at a redshift z = 0.3, and masses M200 = 
0.5, 1 and 2 X 10 h~ Mq. We assume that, because of sub- 
structure in the cluster core and contamination by cluster 
members, the NFW profile is fitted at radii larger than 30 
arcseconds. The inclusion of measurements at smaller radii 
does not change the results significantly, because of the large 
statistical error on small scales. 

Figure ^ shows how the la uncertainty in the determi- 
nation of A/200 (upper panels) and c (lower panels) varies 
with increasing # max , the outermost point included in the 
fit. The uncertainties are determined from 1000 realisations. 
The dotted lines corresponds to the contribution from the in- 
trinsic ellipticities of the background galaxies. The statistical 
uncertainty decreases when the lensing signal is measured 
out to larger distances from the cluster centre. The dashed 
lines denotes the contribution from the large scale structure. 
The combined uncertainty is indicated by the solid lines. 

The results presented in Figure ^ demonstrate that con- 
tribution of distant large scale strucure to the total error 
budget cannot be neglected. This is true, even for clusters 
as massive as M200 = 2 x 10 15 ft -1 Mq. The accuracy in the 
measurement of c improves with increasing mass of the clus- 
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Figure 6. Measurement error in the determination of M200 (upper panels) and c (lower panels) as a function of the maximum fitted 
radius # max for three clusters at a redshift of z = 0.3. The results are obtained by fitting an NFW profile to the simulated tangential 
shear profile, with M200 and c as free parameters. The uncertainties are determined from 1000 realisations. The solid lines indicate the 
total uncertainty, whereas the dashed line denotes the contribution by large scale structure, and the dotted line the contribution from the 
intrinsic cllipticitics of the background galaxies. The results demonstrate that the real uncertainty is a factor ~ 2 larger than expected 
from the noise caused by the intrinsic shapes of the sources alone. Also the accuracy in the parameter estimation does not improve by 
fitting the model to radii larger than ~ 15 arcminutes. 



ter, whereas the uncertainty in M200 increases (although the 
relative accuracy does improve). 

The real uncertainties in M200 and c are typically a fac- 
tor ~ 2 larger than expected from the statistical noise caused 
by the intrinsic shapes of the sources alone. Furthermore, the 
accuracy of the parameter estimation does not improve by 
fitting the model to radii larger than ~ 15 arcminutes. 

It is also interesting to examine how well we constrain 
M200 and c for clusters at various redshifts. For instance, 
Hoekstra (2001) showed that large scale structure seriously 
limits the accuracy of mass determinations for low redshift 
clusters. At higher redshifts, the accuracy is limited by the 
intrinsic shapes of the sources. 

We consider a cluster with M200 = 10 15 /i _1 Mq at red- 
shifts ranging from 0.05 to 0.7. Figure |^ shows that the most 
accurate results are obtained for 9 m « = 15', and we use this 
value to fit the profiles. The la uncertainties in the measure- 
ments of M200 and c as a function of redshift are presented 
in Figure (?[ 

The results show that the optimal redshift for the de- 
termination of A/200 is z ~ 0.3, although for the adopted 
limiting magnitude of R = 26, the uncertainty does not 



increase dramatically for higher redshift clusters. This is be- 
cause for a fixed angular scale, one probes larger physical 
scales at high redshift, which explains why the contribution 
from the intrinsic shapes is almost constant, and the noise 
from large scale structure increases. The noise, however, in- 
creases rapidly for low redshift clusters. 

The situation is different for the measurement of the 
concentration parameter c, which is presented in Figure (7)3. 
The error increases significantly for high redshift clusters, 
and the most accurate measurements are obtained for clus- 
ters with redshifts z — 0.1 — 0.3. 

Therefore studies of intermediate redshift clusters are 
best suited to test the predictions of dissipationless collapse. 
In addition these systems can be studied efficiently using X- 
ray telescopes and new multi-object spectrographs (such as 
IMACS). Comparison of the results obtained through these 
various techniques can provide also new insights in the dy- 
namics of structure formations 
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Figure 7. Measurement error in the determination of Af200 (panel a) and c (panel b) as a function of redshift for a cluster with 
M'200 = W 15 h~ Mq, if an NFW profile is fit to the data with M200 and c as free parameters. The solid lines indicate the total 
uncertainty, whereas the dashed line denotes the contribution by large scale structure, and the dotted line the contribution from the 
intrinsic cllipticities of the background galaxies. The results are for a maximum fitted radius 6> max = 15 arcminutes. The optimal redshift 
for the determination of M200 and c is z ~ 0.3. The uncertainty in A/200 increases slowly with redshift, whereas the uncertainty in c for 
high redshift clusters is large. 



5 CONCLUSIONS 

The comparison of measured mass profiles of cluster of 
galaxies with the outcome of numerical simulations provides 
an important test of our current understanding of structure 
formation. Although mass profiles can be inferred through 
various techniques, weak lensing is of particular interest, be- 
cause it probes the (dark) matter distribution directly. It 
does not rely on assumptions about the dynamical state of 
the cluster or the nature of dark matter. In addition it does 
not require a visible tracer, which is important to study the 
mass profile at large distances from the cluster centre. 

Unfortunately, a simple interpretation of the lensing sig- 
nal is complicated because of the contribution of large scale 
structure along the line of sight. To date, weak lensing stud- 
ies of clusters of galaxies have ignored the effect of distant 
large scale structure. The effect on mass estimates was stud- 
ied previously by Hoekstra (2001). In this paper we followed 
the approach of Hoekstra (2001) to quantify the effect of dis- 
tant large scale structure on the mass profile inferred from 
weak lensing. 

We considered NFW profiles, which provide a good de- 
scription of the outcome of simulations with collisionless 
dark matter. We have created shear profiles, which include 
the noise from the intrinsic shapes of the galaxies and the 
contribution from large scale structure. These profiles were 
fit with NFW profiles, with the mass M200 and the concen- 
tration parameter c as free parameters. 

We found that M200 and c are constrained best for clus- 
ters at intermediate redshifts (z ~ 0.3). For a cluster at 
z = 0.3, optimal results are obtained when the lensing sig- 
nal is measured out to 10 — 15 arcminutes. Measurements 
at larger radii do not improve the accuracy. This is different 



from the situation where the contribution from large scale 
structure is ignored: the measurement error decreases when 
the signal is measured out to larger radii. 

King et al. (2001) studied the effect of substructure on 
the determination of cluster mass profiles, and found that 
the errors in the fitted parameters increase by a few per- 
cent. We found, however, that distant large scale structure 
is a major source of uncertainty in the determination of clus- 
ter mass profiles from weak lensing. The true uncertainties 
in M200 and the concentration parameter c are ~ 2 times 
larger than when large scale structure is not included in the 
error budget. It is also important to note that deeper ob- 
servations do not significantly improve the accuracy of the 
measurements, because the smaller statistical error is coun- 
teracted by an increase in noise due to large scale structure. 



ACKNOWLEDGEMENTS 

I thank Kris Blindert, Ludo van Waerbeke and Howard Yee 
for a careful reading of the manuscript. 



REFERENCES 

Bacon, D., Refregier, A., Clowe, D., & Ellis, R.S. 2001, MNRAS, 
325, 1065 

Bacon, D., Ma ssev. R.. Refregier. A.. & Ellis, R.S., 2002, MNRAS, 



astro-ph/0203134 



submitted, 
Bartelmann, m" 1996, A&A, 313, 697 
Bartelmann, M., & Schneider, P. 2001, Physics Reports, 340, 291 
Bernardeau, F., van Waerbeke, L., & Mellier, Y. 1997, A&A, 322, 
1 

Bonnet, H., Mellier, Y. & Fort, B. 1994, ApJ, 427, L83 



© 0000 RAS, MNRAS 000, 000-000 



Mass profiles from weak lensing 9 



Gioia, I.M. 2000, ApJ, 



Cen, R. 1997, ApJ, 485, 39 
Clowe, D., Luppino, G.A., Kaiser, N 
539, 540 

Clowe, D. & Schneider, P. 2001, A&A, 379, 384 

Clowe, D. & Schneider, P. 2002, A&A, in press, fistro-ph/0208097 

Dubinski, J. & Carlberg, R.G. 1991, ApJ, 378, 496 

Erben, T., van Waerbeke, L., Bertin, E., Mellier, Y., & Schneider, 

P. 2001, A&A, 368, 766 
Fahlman, G., Kaiser, N., Squires, G., & Woods, D. 1994, ApJ, 

437, 56 

Ghigna, S., Moore, B., Governato, F., Lake, G., Quinn, T. Stadcl, 

J. 2000, ApJ, 544, 616 
Gray, M.E., Taylor, A.N., Meisenheimer, K., Dye, S 

Thommes, E. 2002, ApJ, 568, 141 
Hamilton, A.J.S., Kumar, P., Lu, E., 

374, LI 

Franx, M., Kuijkcn, K 



Hoekstra, H 
504, 636 
Hoekstra, H 



Wolf, C. & 
Matthews, A. 1991, ApJ, 
& Squires, G. 1998, ApJ, 



88 



Franx, M., Kuijkcn, K. 2000, ApJ, 532, 
Hoekstra, H. 2001, A&A, 370, 743 
Hoekstra, H., Franx, M., Kuijkcn, K.. 

MNRAS, 333, 911 
Hoekstra, H., Yee, H.K.C., Gladders, M.D. 2002b, NewAR, in 



; van Dokkum, P.G. 2002a, 



press, istro-ph/0205205 
Jain, B., & Seljak, U. 1997, ApJ, 484, 560 
Kaiser, N., & Squires, G. 1993, ApJ, 404, 441 
Kaiser, N., Wilson, G., Luppino, G., Kofman, L.. 



Gioia, I., Mct- 



zger, M., Dahle, H. 1998, ApJ, submitted, istro-ph/9809268 
King, L.J., Schneider, P., & Springel, V. 2001, A&A, 378, 748 
Metzler, C.A., White, M., Norman, M. & Loken, C. 1999, ApJ, 

520, L9 

Moore, B., Quinn, T., Governato, F., Stadel, J., & Lake, G. 1999, 

MNRAS, 310, 1147 
Navarro, J.F., Frenk, C.S., & White, S.D.M. 1995, MNRAS, 275, 

720 

Navarro, J.F., Frenk, C.S., & White, S.D.M. 1997, ApJ, 490, 493 
Peacock, J. A. & Dodds, S.J. 1996, MNRAS, 780, L19 
Schneider, P., van Waerbeke, L., Jain, B., Kruse, G. 1998, MN- 
RAS, 296, 873 

Squires, G., Kaiser, N., Babul, A., Fahlman, G., Woods, D., Neu- 
mann, D.M., & Bohringer, H. 1996, ApJ, 461, 572 
van der Marel, R.P., Magorrian, J., Carlberg, R.G., Yee, H.K.C, 

& Ellingson, E. 2000, AJ, 119, 2038 
Wh ite. M.. van Waer beke. L., & Mackey, J. 2002, ApJ, in press, 



astro-ph/0111490 



Wright, CO. & Brainerd, T.G. 2000, ApJ, 534, 34 



© 0000 RAS, MNRAS 000, 000-000 



